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We construct a real time current-conserving functional renormalization group (RG) scheme on the 
Keldysh contour to study frequency-dependent transport and noise through a quantum dot in the 
local moment regime. We find that the current vertex develops a non-trivial non-local structure in 
time, governed by a new set of RG equations. Solving these RG equations, we compute the complete 
frequency and temperature-dependence of the noise spectrum. For voltages large compared to the 
Kondo temperature, eV S> IcbTk, two sharp anti-resonances are found in the noise spectrum at 
frequencies hw — ±eV , and correspondingly, two peaks in the ac conductance through the dot. 

PACS numbers: 73.63.Kv, 72.15.Qm, 72.70.+m 



Introduction. Not only the current flowing through 
a system, but also its fluctuations (noise) carry crucial 
information on the physics that governs transport 
Zero- frequency noise (shot noise) has been used, e.g., to 
reveal the fractional charge of the quasiparticle excita- 
tions in the fractional quantum Hall liquid Q. The low- 
frequency electrical noise has been extensively studied 
in various systems Q and is by now relatively well un- 
derstood. However, even more information is stored in 
the finite-frequency (FF) current noise: It has been pre- 
dicted that the FF noise is sensitive to the statistics of the 
quasiparticles Q and a crossover between different quan- 
tum statistics can be potentially observed as function of 
frequency, similar to the one observed as a function of 
temperature 0. Moreover, in the quantum regime char- 
acterized by frequencies higher than the applied voltage 
or temperature, the FF noise is a powerful tool to reveal 
the characteristic time scales of the probed system Q as 
well as the dynamics of the excitations or the importance 
of interactions. 

Due to their small size, transport through quantum 
dots (QD) is strongly influenced by the Coulomb block- 
ade. In particular, QDs with an odd number of elec- 
trons behave as artificial magnetic impurities and ex- 
hibit the Kondo effect Q, a paradigmatic many-body 
phenomenon corresponding to the screening of the spin 
of the quantum dot by the conduction electrons of the 
leads at temperatures T below the Kondo temperature, 
Tff. QDs thus provide an ideal test ground to study 
non-equilibrium transport in the presence of strong in- 
teractions. While the conductance of a QD in the Kondo 
regime is well understood by now Q , much less is known 
about current-fluctuations. Though these are promising 
quantities to characterize the out of equilibrium Kondo 
effect, nevertheless, most experiments focused so far on 
measurements of the average current , and even 

results on low-frequency noise measurements have only 



appeared recently [12| . While the effect of ac voltage on 
the non-equilibrium Kondo effect has been studied ex- 
perimentally relatively longtime ago [13j ]. no FF noise 
measurements have been reported so far to our knowl- 
edge. Theoretically, most studies focused on shot noise: 
a non-monotonous bias-dependence of the shot noise with 
a maximum at eV ~ fcsTif has been found at T <C Tk 
[h3 |. and a universal ratio 5/3e between the shot noise 
and the backscattering current at T = has been pre- 
dicted for the SU{2) Kondo effect [l5|. 

The purpose of this work is to provide a general 
analysis of the finite-frequency current noise through a 
quantum dot in the local moment regime. To achieve 
this, we construct a real time functional renormaliza- 
tion group (FRG) scheme on the Keldysh contour to 
study frequency-dependent transport and noise through 
a Kondo quantum dot. Our formalism reproduces the 
scaling equations of Rosch et al. [l6[ for the vertex func- 
tion. However, we find that the current vertex also de- 
velops a non-trivial non-local structure in time, governed 
by a new set of RG equations. Such structure of the 
current vertex turns out to be unavoidable to guaran- 
tee current conservations and is necessary to calculate 
the finite frequency current noise in a controlled manner. 
Solving this set of RG equations, we compute the com- 
plete frequency and temperature dependent noise spec- 
trum through the dot. Our approach is valid at any fre- 
quency lu, voltage V and temperature T provided that 
maxjeV, fc^T} > HbTk- For frequencies huj IcbTk, 
we find sharp anti-resonances in the voltage dependence 
of the noise spectrum at eV — fiio, which gradually disap- 
pear with increasing temperature. The absorption noise 
is also found to exhibit strong anomalies at fiw = eV, and 
we find similar anomalies in the non-equilibrium ac con- 
ductance, too, where a split non-equilibrium Kondo reso- 
nance is observed. Precursors of the noise anomaly have 
been found in the zero-temperature symmetrized noise at 
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finite frequency, as first computed at the Toulouse point 
of the Kondo model [l7|, and later confirmed by a non- 
equilibrium one- loop perturbative calculation [18[ . How- 
ever, logarithmic singularities are completely absent at 
the rather special Toulouse point 17J , while the method 
of Ref. [lH was not accurate enough to capture fine de- 
tails of the anomaly. 

Model. In this paper, we focus our attention to the lo- 
cal moment regime of the quantum dot, where we can de- 
scribe the electron on the dot as a spin 5 = 1/2 moment 
which couples to electrons in the left and right electrodes 
through an exchange interaction Q, 
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Here a stands for the three Pauli matrices, the fields 
V-W = / c Q(T (£)e - ^' d£ destroy electrons of spin a in 
leads a <E {L,R}, with 1/a a high energy cut-off (2f| . 
The dynamics of i\) aa is governed by the non-interacting 
Hamiltonian, H Q = X) Q>(T /d£ (f + Ma) 4 CT (0 c aa(6> 
with fi a = eV a , the chemical potential shift of lead a. 

To describe the spin using standard field-theoretical 
methods, we make use of Abrikosov's pseudo-fermion rep- 
resentation [l9| : we introduce a fermion operator /j for 
each spin component, s = ±1/2, and represent the spin 
operator as S" — > J2 S s' \ fs a l with the additional 
constraint, J2 S fjf s = 1. 

We then employ a path integral formalism on the 
Keldysh contour. In this approach each fermionic field 
is replaced by two time-dependent Grassmann fields liv- 
ing on the upper and lower Keldysh contour (n = 1,2), 
respectively, and the dynamics is determined by the 
Keldysh action, S = S\ ca d + <S sp in + Sint • The parts <Si ea d 
and iS sp in describe the conduction electrons and the spin 
in the absence of interaction. They are quadratic in the 
fields, and determine the non-interacting Green's func- 
tions m. 

The interaction part of the action, iSj nt , is diagonal in 
the Keldysh indices and is initially local in time. How- 
ever, elimination of high energy degrees of freedom in 
course of the RG procedure generates retardation effects, 
and the interaction becomes non-local. We find that, 
with a good approximation, it can be expressed as 

Sint = EE Sk ^ / dil d ^ 2 9a,p(h - t 2 ) 

f K (T 12 )ar(T 12 ) ■ C(ti)o^(i 2 ) , (2) 

where T\ 2 = (t\ + t 2 )/2 and s K — ±1 for the upper and 
lower Keldysh contours, respectively. The initial (bare) 
coupling function g a p{t) is local in time, and is given by 

ffoflW = 3 a P ^(*)* The justification for this structure, 
Eq. (f2|) , is straightforward: The spin evolves very slowly, 
and its time evolution can be very well approximated by 
the one in the absence of interactions at electronic time 



scales. However, conduction electrons have fast dynam- 
ics, and their retardation effects become important as one 
approaches smaller energy scales. 

Functional RG. We construct the RG equations by 
expanding the action in S- ln t and rescaling the cutoff pa- 
rameter a — > a' . An intcgro-differential equation is ob- 
tained for the functions g a p(t), which becomes simple in 
Fourier space, 



dl 



= gM q(w,a) g(w) 



(3) 



Here I = ln(a/ao) is the scaling variable, ao is the ini- 
tial value of the cut-off time, and we introduced the ma- 
trix notation, g a p — > g. The matrix q(w, a) is a cut-off 
function, which depends somewhat on the precise cut-off 
scheme, but for practical purposes is well-approximated 
by the function q a g(uj, a) w 8 a p 0(i — |a; — fjb a \) at 
T = temperature [l(|. The scaling equation, Eq. 
is identical to the one obtained in a more heuristic way 
in Ref. [lij . however, in our real time functional RG for- 
malism the derivation is rather straightforward and sim- 
ple [2l[. We remark that the usual poor man's RG proce- 
dure can be recovered by dropping the time-dependence 
of g a p, and replacing the generated non-local couplings 
by local ones, g a p(t) ~ * S(t) J dt g a p(t), which corre- 
sponds to setting us — >• in Eq. (j3]). 

Our primary purpose is to compute current-current 
correlation functions. To do that, we first define the left 
and right current operators from the equation of motion, 

giving i L (t) = -i R (t) = E a/J §«&s(i) • ft a (t)°Mt), 

with the current vertex matrices defined as 



v L = -v R = 
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In the path integral language, it is useful to introduce 
a corresponding generating functional, 



ZK(t)] = (e 



(5) 



from which the current-current correlation functions can 
be generated by functional differentiation with respect 
to h*(t). A systematic investigation of the leading loga- 
rithmic diagram series shows that the expression of the 
current field, I K (t), necessarily becomes non-local in time 
under the RG procedure, and acquires a form, 



7 x(*) = I EE / di 1 dt 2 V^(ti-i,t-ta,o) 

K Oip 
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(6) 



The physical motivation of the double time-structure is 
simple: in the renormalized theory it is not enough to 
know the times electrons enter and leave the dot {t\, 2 ), 
but the time of the current measurement, t, must also be 
kept track of. 
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FIG. 1: (color online) Voltage-dependence of the symmetrical 
noise, as computed by FRG for ftw = 30 UbTk- 



It is relatively straightforward to derive the scaling 
equations from the perturbative expansion of Eq. ([5]). 
and we obtain 



dV L ((J!,M2) 

dl 



= V L (wi,w 2 ) q(wa,o) g(w 2 ) 

+ g(wi) q(wi,o) V L (wi,tj2) ■ (7) 



This equation needs be solved parallel to the scal- 
ing equation, Eq. ((3|) with the boundary condition, 
V L/R (T 1 ,T 2 ,a ) = 5(n) <5(t 2 ) v L / fl . Though the renor- 
malized couplings gz,(u;) drive the scaling of the cur- 
rent vertexes, Vi(wi,a; 2 ,ao), there seems to be no sim- 
ple connection between these too. In other words, it is 
unavoidable to introduce the renormalized current ver- 
texes within the functional RG scheme to compute time- 
dependent current correlations. The above extension 
seems to be also necessary to guarantee current conser- 
vation: Eq. is linear in V L , and therefore the con- 
dition I1(t) + Ifjit) = is automatically satisfied, for 
any value of the cut-off, a. On the other hand, we could 
not find any way to generate a current field from just 
the renormalized action, Eq. ((2J, such that it respects 
current conservation. 

We solved Eqs. ([3]) and (O numerically to obtain g 
and W L / R . Both display singularities at frequencies 
huj = ±eV/2 [2(J. With the couplings g(wi) and the cur- 
rent vertexes \ rL ^ R (u>i, w 2 ) in hand, we then proceeded to 
compute the noise through the device, by doing perturba- 
tion theory with the renormalized action. For the Fourier 
transform of the absorption and emission noise compo- 
nents S> L (t) = (I L (t)I L (0)), and S< L (t) = (I L (0)I L (t)), 
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we obtain S^ L (u>) 



s ll(-u), with 



£iM = y S(S + l) J ^Tr{V 

G > {u + )V L (u + ,Q-)G<(G>-)} 



UJ+) 



(8) 



Here Q± = ui ± and the bigger and lesser Green's 

functions G^ < (w) = ±i 2ir S a p — Ma))- 

Results. The symmetrized noise spectrum, Sll{^>) = 
\[Sl L {uj) + S^ L {uj)] is plotted in Fig.fflfor huj = 30 k B T K 




FIG. 2: (color online) Temperature-dependence of the sym- 
metrical noise, as computed by FRG. Inset: Comparison with 
third order renormalized perturbation theory (RPT) using a 
reduced cut-off, a — h/WeV and the corresponding couplings 
jLn(ci), and effective second order perturbation theory with 
renormalized couplings jii?(a = h/eV), but still using the 
original bandwidth. None of these methods are able to get 
close to the FRG results. 



as a function of voltage, V. Clearly, the noise spec- 
trum shows rather strong features at the bias voltages, 
eV Huj. The appearing dips are clear fingerprints of the 
non-equilibrium Kondo effect, and they gradually vanish 
as we increase the temperature, T [26|. The increase in 
the noise with decreasing V at low voltages can be under- 
stood as being the consequence of increasing spin relax- 
ation time, resulting in an increased Kondo conductance 
through the dot. At high voltages, on the other hand, 
increased photon emission is mainly responsible for the 
increasing noise, thereby giving the V-shaped pattern. 

Similar features appear in the frequency-dependent 
symmetrized noise (see Fig. [2]). It is instructive to com- 
pare the FRG results with perturbation theory, giving 



si L {t) 



s 2 jliifll 



cos(eVY) j- 



1 



+ 2(j LL +j RR ) 



(*■ 

ln(l + i t/a) 



iaf 



•}■ (9) 



t(t-2i a) 

The curves in the inset of Fig. [5] are obtained by tak- 
ing the Fourier transform of this expression. While 
the perturbative result also exhibits singular features at 
huj = ±eV, however, it does not reproduce the precise 
shape of the anomaly, even if we use renormalized pa- 
rameters, j a p —t j a p(eV) obtained by solving the usual 
leading logarithmic scaling equations. 

Experimentally, it may be more convenient to measure 
separately the emission or absorption noise components 
at a fixed finite frequency, S P /„ (to > 0) = S < (±lj), as 
function of the bias voltage |22|. As shown in Fig. [3l at 
T = the emission noise vanishes at voltages V < hcu/e 
due to energy conservation, and has an abrupt logarith- 
mic singularity at V = fko/e at temperatures T <C Tk, 
which is gradually smeared out for fc^T > eV. 
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quantities which are within experimental reach 22 1. 
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FIG. 3: (color online) Voltage and temperature-dependence 
of the emission noise as computed through FRG for fku — 
5k B T K . 



Another quantity which may be easier to access exper- 
imentally is the non-equilibrium finite frequency linear 
conductance, defined as the current response of the sys- 
tem to an external time-dependent variation of one of the 
lead potentials. According to a formula of Safi [23J, this 
can be expressed as 



Re G ll (lu,V) = j-(S^ L H 



(10) 



Notice that Gll{0, V) is just the usual non-equilibrium 
differential conductance, Gll(0,V) = &v l {Il)^ while 
Gll(u, V — > 0) corresponds to the usual equilibrium ac 
conductance [24j . Re Gll(u, V) is an even function of ui 
and exhibits two peaks at hui = ±eV, associated with the 
non-equilibrium Kondo effect [l4| (see Fig. @| . This con- 
firms that the dips in the noise are related to the splitting 
of the Kondo resonance at a finite bias. 

Summary. We have developed a real-time functional 
rcnormalization group approach. We have shown that 
the current vertex becomes non-local in time under renor- 
malization which turns out to be necessary in order to en- 
sure non-equilibrium current conservation. Within this 
formalism, we have been able to calculate the voltage 
and temperature dependence of the current noise at fi- 
nite frequency and the non-equilibrium ac conductance, 
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FIG. 4: (color online) The ac conductance for ui > 0. 
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